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We re-examine the band structure of the stripe charge ordered state of a-(BEDT-TTF)2l3 under 
pressure by using an extended Hubbard model within the Hartree mean-field theory. By increasing 
pressure, we find a topological transition from a conventional insulator with a single-minimum in 
the dispersion relation at the M-point in the Brillouin zone, towards a new phase which exhibits 
a double-minimum. This transition is characterized by the appearance of a pair of Dirac electrons 
with a finite mass. Using the Luttinger-Kohn representation at the M-point, it is shown that such a 
variation of the band structure can be described by an effective 2x2 low energy Hamiltonian with 
a single driving parameter. The topological nature of this transition is confirmed by the calculation 
of the Berry curvature which vanishes in the conventional phase and has a double peak structure 
with opposite signs in the new phase. We compare the structure of this transition with a simpler 
situation which occurs in two-component systems, like boron-nitride. 
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I. INTRODUCTION 

Two-dimensional molecular conductor q;-(BEDT- 
TTF)2l3i has brought much interest by the variety of 
electronic states, such as an insulating stripe Charge Or- 
dered (CO) stato^, a superconducting state in the pres- 
ence of charge ordering, and a Zero Gap State (ZGS) 
with a masslcss Dirac spectrum^. 

The stripe CO state^i^i^, which was suggested to ex- 
plain an insulating phase below 135K at ambient pres- 
sure, was confirmed by NMR experiment^. The super- 
conducting state, found under uniaxial pressure along 
stacking axis (a-axis)^, was investigated theoretically by 
using an extended Hubbard model^. A narrow gap state 
was suggested by Kajita to explain anomalous increase 
of Hall coefficient at high pressuresi^. The existence 
of massless Dirac electrons was predicted theoretically^i 
based on a tight-binding calculation using the transfer 
energies of ref. [l^- Thus it has been revealed that 
this narrow gap state is indeed a zero gap state, as also 
confirmed by a first principle calculatio n^^d^ . In this 
state, the energy spectrum near the Fermi energy con- 
sists of two cones described by a tilted Weyl equation 
for massless Dirac electron o^^'^^ . This has been con- 
firmed by a comparison between the theoretical and ex- 
perimental results for the temperature dependence of the 
Hall coefficientiii^ and the angular dependence of the 
magnetoresistancei^ii^. 

The two tilted Dirac cones were predicted to merge and 
disappear at the F point under extremely high pressure^^. 
This merging transition was also studied in the context of 
deformed graphene with a variation of transfer integrals, 
and was described using a generalized two-component 
Hamiltonian for Dirac electrons. This effective Hamil- 
tonian describes the merging of two Dirac cones and the 
opening of a gap at the transition^SiSii^. 

In the present paper, we show that another type of 
transition may also occur in the CO state of q;-(BEDT- 



TTF)2l3. In this state, the energy spectrum exhibits a 
gap between the conduction and valence bands at the M- 
point, with a single minimum. By increasing pressure, it 
is demonstrated that new electronic phases emerge in the 
CO state with a double- minimum structure in the vicin- 
ity of the M-point. This double- minimum corresponds 
to the emergence of a pair of massive Dirac electrons, 
whose study is the main goal of this paper. This qual- 
itative change of the band structure is described by an 
effective Hamiltonian with a single driving parameter. 



The paper is organized as follows. In section 2, a tight 
binding model for a-(BEDT-TTF)2l3 is described where 
the repulsive interactions between molecules are treated 
within Hartree mean-field theory. When varying pres- 
sure and the intersitc repulsive interaction, a new phase 
diagram is obtained, and is described in section 3, where 
new phases characterized by a Dirac electron pair are ob- 
tained in the CO state. In section 4, using the Luttinger- 
Kohn representation at the M-point of the Brillouin zone, 
we construct a 2 x 2 effective Hamiltonian to describe the 
low energy band structure near this point, and the emer- 
gence of a pair of Dirac points. The general structure of 
this effective Hamiltonian also describes the emergence of 
Dirac points in a simple toy-model related to the physics 
of boron nitride (BN). In section 5, we show that this 
Dirac pair is revealed by two sharp peaks of opposite 
sign in the Berry curvature^^. Then each Dirac point is 
characterized by an appropriate Berry phase calculated 
by using a method already applied to the physics of Dirac 
points in boron nitride^i. The last section is devoted to 
summary and discussion. 
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FIG. 1: The model describing the electronic system of a- 
(BEDT-TTF)2l3i^'2S^. The unit cell consists of fom BEDT- 
TTF molecules A, A', B and C with seven transfer energies. 
The nearest neighbor repulsive interactions are given by Va 
and Vb- The a- and &-axis in the conventional notation corre- 
spond to the y- and x-axis in the present paper. 



II. TIGHT-BINDING MODEL FOR 
a-(BEDT-TTF)2l3 

A. The Hamiltonian 

The model used to describe the two-dimensional 
electronic system in a-(BEDT-TTF)2l3 is shown in 
Fig. B^i^s^. The unit cell consists of four BEDT-TTF 
molecules on sites A, A', B and C. The sites A, B and 
C are inequivalent, while A is equivalent to A' so that 
inversion symmetry is preserved. There are six electrons 
for the four molecules in a unit cell, i.e., the band is 
3/4-filled. On the basis of the HOMO orbitals of these 
sites^i^, these electrons are described by the extended 
Hubbard model with several transfer energies, the on- 
site repulsive interaction U and the anisotropic nearest- 
neighbor repulsive interaction Vap : 

H = Yl (*»";J/3 O-La^jPa + Inc.) 

{ia:jf3),a 

ia 

+ Yl ^"/^ alaa^]l3a'"'3fl'^'"'ic<y^ (2-1) 

{ia:jP).(7,(T' 

where i,j denote site indices of a given unit cell, and 
a,/3(= A, A', B and C) are indices of BEDT-TTF sites 
in the unit cell. Hereafter, the energies are given in eV. 
In the first term, a]^^ denotes a creation operator with 
spin cr(=t, i) and tia-jp is the transfer energy between 
the (i, a) site and the (j, /3) site. 

The transfer energies as a function of a uniaxial pres- 
sure (Pa) along the a-axis are estimated from an extrap- 



txiPa)^txmi + KxPa). 



(2.2) 



The transfer energies tx and the coefficients Kx are ob- 
tained from the numerical data at Pa = Okbar— i^ and 
at Pa = 2kbari^, such that iai(O) = -0.028, <q2(0) = 
0.048, ia3(0) = -0.020, i6i(0) = 0.123,42(0) = 
0.140,^3(0) = -0.062,44(0) = -0.025, [eV] and Kai = 
0.089, Ka2 = 0.167, Ka3 = -0.025, A'bi = 0,Kb2 -= 
0.011, A'h3 = 0.032, Km = [eV/kbar]. We use the pa- 
rameter U = 0.4 and the Vap take two different values, 
Va = 0.17 ~ 0.18 eV along the stacking direction, and 
Vb = 0.05 eV along the perpendicular direction^. With 
this choice of parameters, we obtain a pressure depen- 
dence of the electronic spectrum consistent with exper- 
imental results^ii^. Moreover the results of the mean- 
field theory coincide with the experimental results for 
the charge disproportionation in the ZGS and the charge 
ordering in the CO phase^ii^. By "charge dispropor- 
tionation" , we mean that the A, B and C sites in a unit 
cell are inequivalent but the inversion symmetry between 
A and A' remains. By "charge ordered" state, we mean 
that the inversion symmetry between A and A' is broken. 
Throughout the paper, h and the lattice constant a are 
taken as unity. 



B. Hartree mean- field theory 

As in previous work o^i^^ , we restrict ourselves to a 
Hartree mean field theory which implies that the mean 
field Hamiltonian is diagonal in spin space. In addition 
we only consider mean field solutions that do not break 
the underlying Bravais lattice symmetry; this further im- 
plies that for a given spin a the mean field Hamiltonian 
Her (k) is a 4 X 4 matrix in the Bloch-basis Ok^a . In second 
quantized form, the mean-field Hamiltonian reads : 



H 



MF — 



BZ 



(2^ 



E^-w, 



aP 



(2.3) 
(2.4) 



ea/3cr(k) = Iac,5ap+eap{\<i), (2.5) 



P'a' 



ea/3(k) = yjape'^'^i 



(2.7) 



,t 



where (uaa) = {O'laa^iacr) is the mean field local density 
of spin a for the molecular state a in unit cell i; {uaa) 
is assumed to be independent of the unit cell i. 5 de- 
notes the vectors connecting nearest neighbors sites. The 
spin dependent site potential laa represents the Hartree 
mean field which comes from the on-site Hubbard U and 
nearest- neighbor Coulomb interactions Va,Vb. 



Note that ofF-diagonal elements eQ,^(k) are indepen- 
dent on spin index a and on interaction strength but 
depend on k and Pa- Moreover these off-diagonal ele- 
ments ea^(k) have the time reversal symmetry e*^(k) = 
£(1/3 (—k). In this representation, only the self-consistent 
diagonal elements may explicitly exhibit the breaking of 
time reversal and/or inversion symmetries. 

To obtain the mean-field phase diagram, the Hamil- 
tonian is diagonalized numerically for a given k in each 
spin subspace, according to 



^ea/3<T(k) d/37<T(k) = ^^.^(k) da^<^(k), (2.8) 

where S,-ya are the eigenenergies ordered such that, 
Cia(k) > 6.(k) > e3a(k) > e4.(k) (7 = 1,2,3,4 is the 
band index), and da^ai}^) are the corresponding eigen- 
vectors. In terms of cq. (|2.8p . the average number {uaa) 
of electrons with spin cr on a type of site is expressed as 



0.2 



dk 



BZ 



(27r) 



2 ^'^a^<yO^) dajc,{k)nF{£.ja{^)), 



7 = 1 



(2.9) 



where npi^-yaCk)) = l/(exp [(^^^(k) - ^i)/T] + 1) is the 
Fermi factor at temperature T (fcs — 1) and fi is the 
chemical potential determined from the condition of a 
3/4 filled system: 



3 

4 



^E("-) - ^X^ (^E"^(^7.w). 



(2.10) 



Equation p.9p constitutes the self-consistent relation for 
the mean field quantities (riaa-)- A high accuracy in ful- 
filling these constraints requires a very fine mesh of the 
reciprocal space which constitutes the main numerical 
difficulty. When convergence is achieved, the mean field 
energy per unit cell of a given state is then calculated as 



E 



MF 



/ T^^^Fii'T'^i^M-f'^i^) 



- ^Ua{na-l){nal) - ^ Va/3 (riacr) (»^0<t' ) • 

(2.11) 
The ground state is obtained from minimization of -Emf- 
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FIG. 2: The phase diagram on the uniaxial pressure along a- 
axis {Pa [kbar]) and the repulsive interaction between nearest- 
neighbor sites along a-axis (Va [eV]) where U = 0.4 eV, and 14 
= 0.05 eV. The CO and COM denote insulating and metallic 
states respectively. In addition to the phase (I) of the previous 
work^, there exists a new phase (II), which is characterized 
by a double minimum in the up spin band. 





III. PHASE DIAGRAM 

The (Pa-Va) phase diagram obtained from the self- 
consistent Hartree approximation described above is pre- 
sented in Fig. [5] Pa, is the uniaxial pressure along the a- 
axis and Vq is the repulsive interaction between nearest- 
neighbor sites along the a-axis. This figure exhibits three 



FIG. 3: Schematic behavior of the energy spectrum close 
to the Fermi energy (horizontal line), in the different sub- 
phases of the CO phase : (a)CO(I), (b)CO(II), (c)COM(I), 
(d) COM(II). The red and blue bands correspond respectively 
to t and ^ spins. The center of the horizontal line corresponds 
to the M point (kji/ = (7r,±7r)). 



transition lines. Two of them (the continuous and dashed 
hues) were already found in previous works^ and the third 
one (dotted line) is a novel transition that constitutes the 
main object of this work. The schematic band spectrum 
close to the Fermi surface is shown in Fig. [3] and is dis- 
cussed below. Before describing this last transition in 
more details, we remind the qualitative picture associ- 
ated to the two other transitions. 

In Fig. [21 the solid line marks a charge ordering transi- 
tion resulting from the simultaneous breaking of time re- 
versal and inversion symmetries. It separates the charge 
ordered metallic state (COM) from the Zero Gap state 
(ZGS). In the ZGS phase, energy bands are spin degen- 
erate and inversion symmetry is preserved such that con- 
duction band ^i^. (k) and valence band ^2cr (k) touch each 
other at two Fermi points (Dirac points) k-t . These Dirac 
points move in the Brillouin zone when varying param- 
eters Pa and Va- Around each point kj-, the dispersion 
relation is linear. In addition there is a large anisotropy 
in the Fermi velocities (a factor ^ 10 between the high- 
est and lowest velocity value a^^i^^ ), this appears as a tilt 
in the Dirac cones^^. Coming from this high pressure 
ZGS phase and traversing the continuous line, the in- 
version symmetry is spontaneously broken by the elec- 
tronic interactions. As consequence, for a given spin a, 
a gap opens between bands ^i^{]<.) and ^20- (k), leading a 
priori to an insulating phase. However the time rever- 
sal symmetry is also spontaneously broken by the inter- 
actions so that the degeneracy between t and I bands 
is now lifted. Therefore the simultaneous breaking of 
time reversal and inversion symmetries results in a semi- 
metallic phase (COM) with band overlap leading to small 
electrons and holes pockets of opposite spin orientations 
(Fig.©. 

In striking contrast with the continuous line, the 
dashed line marks a metal-insulator transition from a 
charge ordered metallic (COM) phase to a charge or- 
dered insulator (CO) without breaking of any symmetry. 
In traversing this transition line, the dispersion relations 
of the four energy bands £,ia(i<-) and C2o-(k) stay simi- 
lar but their relative positions to the Fermi level vary in 
such a way that, in the CO phase, the Fermi level falls 
in a true charge gap that separates a valence band and 
a conduction band of equal spin orientation (Fig. |3]). We 
stress that, at the metal-insulator transition, the elec- 
tron density (riaa-) and therefore the resulting mean field 
potential I^a exhibit a cusp on respective sites. 

In this work, by a more detailed analysis of the COM 
and CO phases, we find a new topological transition (dot- 
ted line in Fig. [2]) that further splits each of the COM and 
CO phase into two phases: COM(I,II) and CO(I,II). The 
electronic spectrum in each phase is represented schemat- 
ically in Fig. [3] . The band structure of both t band 
and I band in CO (II) state is explained in Appendix A. 
This transition concerns a modification in the two energy 
bands close to the Fermi energy. They correspond to a 
given value of the spin, that we choose to denote by t- 
The two other bands {],) are not concerned by this tran- 



sition. As illustrated in Fig.|4l the transition from CO (I) 
(Fig.|3]a) to CO(II) (Fig.|3]b) is characterized by a change 
in the form of the dispersion relation of the valence and 
conduction bands. In the CO (I) phase, there is a single 
minimum of charge gap whose position in k space stays at 
the M-point k^v/ = (tt, ±7r), independently of the param- 
eters Pa,Va. Around this point kyv/, valence and conduc- 
tion bands disperse quadratically. In the CO (II) phase, 
the single charge gap separates in two points at symmet- 
rical positions (For example, k± = ±(0.95, — 0.71)7r for 
Pa = 5.4 kbar and T4 = 0.18 eV) from the M-point, and 
k-t move continuously with parameters Pa,Va. There is 
now a double-well structure in the dispersion relation. 
The aim of this work is to describe this topological tran- 
sition in the framework of a universal Luttinger-Kohn 
Hamiltonian. 

As a last remark concerning the phase diagram of 
Fig. [21 we note an intersection of the dashed line and the 
dotted line, at the critical point of coordinates {P*, V*) 
=(4.75 kbar, 0.175 eV). Such crossing implies the prox- 
imity of several minima of the total energy in its neigh- 
borhood. Therefore in this region of the phase diagram 
one may need to go beyond Hartree mean-field approx- 
imation. This might be an important issue since this 
critical point P* , V* is in a region of parameters value 
that is reachable experimentally. 



IV. TOPOLOGICAL TRANSITION AROUND 
THE M-POINT 

A. Low energy 2x2 effective Hamiltonian at the 
M-point 

In order to analyze the band structure in the CO state, 
we expand the mean-field Hamiltonian up to the second 
order in momentum around the M-point, writing k — 
ki/ + Q- The Hamiltonian is diagonal in spin space and 
it has the 4x4 structure, for each spin direction a : 



HaPaC^) — HaPaO^-Al) 



1 



+g^a^ff„/3,(k)|k,, + -q^'q''^^^,Ha,p^{k)\^,,, 

H^f3a{kM)+SH^f3a{<l)- (4.1) 

We wish to restrict ourselves to the two bands which are 
close to the Fermi level (ct =t) ■ This can be done by using 
the Luttinger-Kohn representation^ at the M-point, and 
we obtain an effective 2x2 Hamiltonian with the matrix 
elements {h',h'' — 1,2) up to O(q^) ;^- 



H, 






1 



>U Six" Si/' Si>" 

(4.2) 

where ^° and |j/) are respectively the eigenvalues and the 
eigenvectors of HCknj)- Thus the effective Hamiltonian 



0.02 




-0.02 



FIG. 4: The q-dependence of the conduction and valence 
bands for t spin where q = k — Icm, and the center denotes 
the M point (kj\/ — (vr, ±7r)). a) Spectrum of the CO(I) phase 
with a single minimum at the M point, b) Spectrum of the 
CO(II) phase with a double minimum around the M point. 
Here the parameters are Va = 0.18 eV and Vb = 0.05 eV, and 
the pressures are respectively Pa = 4.5 kbar and Pa — 5.4 
kbar, for the CO (I) and CO (11) phases. 



is rewritten as 

iJ^^'(k) = H"^kM) + tq. ■ V(kM) + Y. W^J{kM)q^qJ, 

ij 

(4.3) 
where, by construction, H^^Ck^f) is diagonal. Like the 
original Hamiltonian, the Luttinger-Kohn (L-K) Hamil- 
tonian has the symmetry H^^ (k) = iJ-^-^(— k). There- 



fore, near the M-point as well as in the vicinity of 
the other so-called time-reversal points of position G/2 
where G is a reciprocal lattice vector, the energy bands 
^a (k) have peculiar properties : for any band ^a (k) , the 
particular values ^q,(G/2) are either a local extremum or 
a saddle point. As a consequence, the gap between two 
consecutive bands at these symmetry points is also either 
a local extremum or a saddle point. This picture explains 
the topological transition from the CO(I) phase where 
the gap at M is a local minimum to the CO (II) phase 
where the gap at M is a local saddle point^. To be more 
precise, let us focus on the structure of the Luttinger- 
Kohn Hamiltonian (|4.3p . A strong consequence of the 
symmetry H^^*{'k) = 7J^^(-k) is that at the M-point, 
like at the other symmetry points, Wij (k^f ) are real sym- 
metric matrices and T^(kM) are real antisymmetric ma- 
trices. We deduce that the minimal form of the 2x2 
matrix H^^Ck^j + q) to order qiQj can be recast as: 

i7^^(kM+q) = /o(q)^o + /i(q)ai+/2(q)^2+/3(q)^3 , 

(4.4) 
where we have used a representation in terms of Pauli 
matrices at and 



/o(q) = A^ + X^'^o' (i^ij ' 

ij 

/i(q) = X]^i^ *^j ' 

ij 

/2(q) = v-q, 

/3(q) = A + ^Wg^ qtqj . 



(4.5) 

(4.6) 

(4.7) 
(4.8) 



with A > 0. The eigenenergies near the M-point are 
given by 



^±(q) = /o(q) ± VhW+hW+hin) 



(4.9) 



The local band structure around M has the symmetry 
^±(q) ~ C±(^q)- Moreover these expressions also imply 
that dq-^±{q)\q=o = so that, as anticipated, M is ei- 
ther a local extremum or a saddle point of the dispersion 
relation of each band. The same properties are also valid 
for the gap A(q) = ^(C+lq) ^ '?-(q)) separating the two 
bands. 



B. CO(I)-CO(II) transition: phase boundary 

Within the framework of this Luttinger-Kohn repre- 
sentation, we now determine the condition for a transi- 
tion from a minimum to a saddle point at the M-point. 
For this purpose, it is useful to define the 2x2 stability 
(Hessian) matrix Sm. We have 



A(q) 



C+(q)-e-(q) 



A- 



l^iqx qy)SM 



(4.10) 



where 



To lowest order in q, the gap function can be expanded 
as : 



Sm = 



with ■ 



O'A(q) a"A(q) 



"9r 



a„., d„. 



3"A(q) a"A(q) 



M 



Wo 



(4.11) 



-TTir- ■ The determinant of this matrix: 



2A 



det^M = w^wf^ 



/ — xy\2 



(4.12) 



governs the stabihty of the M-point. There is an ex- 
tremum when det5'j\/ > and a saddle point when 
detS'M < 0. Therefore the transition line is given by 
the condition det^M = 0, that is 



A(q) ^ A 



2A 



)9^. 



(4.17) 



In the CO(II) phase, the function A(q) has a saddle 
point. This corresponds to the case where detSM < or 
equivalent ly Wg^Wg^ — (Wg**)^ < 0. There is an angular 
region where A(q) < A. This happens in an interval 

2 

0min < ^ < Smax sucli that (wg + ^) < 0. The two an- 
gles 0min, Omax, choscu in the range [0, tt], are determined 
by the condition (uig + 5^) = 0, which gives 



tan 6*, 



-Wg^ ± V— detS'M 



max .min 



(4.18) 




-(-3^ + lS^ 

(4.13) 
We emphasize that the stability matrix Sm is totally 
independent of the terms wl^ . Nevertheless, these terms 
are important to determine the position of the Dirac 
points and the size of the gap at the Dirac points in 
the CO(II) phase. They also play a crucial role in deter- 
mining the topological properties of the band structure 
(Berry curvature). Therefore, among the 9 initial param- 
eters of Eqs. (|4.6p . (|4.7p and (|4.8|) , only 6 are pertinent 
to determine the transition line. They enter the combi- 
nation detS'M which is the single driving parameter for 
the CO (I) - CO (II) transition. 



C. Minimal form of the effective Hamiltonian for 
the merging transition 

In order to analyse the local structure of the Hamil- 
tonian near a time-reversal point, it is convenient to 
parametrize it using polar coordinates of the wave vector 
q. We write : 



qx = q cos f 



qy = q sm f 



(4.14) 



so that the Hamiltonian describing the CO(I)-CO(II) 
transition can be written in the form (since the compo- 
nent ctq pl^'YS no role in our discussion, we define the ef- 
fective Hamiltonian as h{q) = H^^"^ (kM + (l) — /o(q)fo) ■ 



Mq) = 



A + w^q^ —ivoq + wfq"^ 
ivgq + wfq"^ —A — w^q"^ 



(4.15) 



where 



wf = vof SIT? 6 + wf'^ cos^ 6 
Wg = lUg^ sin^ 9 + wf ^ cos^ 9 
vg = Vx COS 9 + Vy sin 9. 



(4.16) 



The gap function can be rewritten in terms of these an- 
gles as 



A(q) = A + w sin{9 - 9max) sin(6' - 9„itn)q^ , 

^[cos6'_-cos2(6'-6'+)]g2 (4.19) 



A- 



where w can be rewritten conveniently as 

,1/2 



b 



'M 



4detS' 



M\ 



WM = \wT'Vwf\ 



The angles 9^ and 9- are defined as 



1 







(4.20) 



(4.21) 



They represents respectively the most unstable direction 
and the angular aperture of the region A(q) < A. Fig. 
[5ja) shows the gap function in the CO(II) phase at Pa = 
5.4 kbar. Their expression is 



TT 1 / 2I(J?^ 

- + - arctan ( _,, _ _yy 



Wo 



Wo 



arctan 



2V-detS'. 



M 



WM 



(4.22) 



Near the merging transition, when Pa > Pm , &+ -^ 6m 
where the angle 9m given by 



n[ 



TT — arctan 



TTivy 



(4.23) 



defines the direction of emergence of the two Dirac points. 
Moreover if — > wm — \w'f^ + w^ \ and the angular aper- 
ture 9- vanishes as 9^ ~ 2 V— detS'M /wAf. 

Near the transition, we define a new set of coordinates 
along the direction of emergence and the perpendicular 
direction. Writing q'^ = qcos{9 — 9m) and q'y = qiiin{9 — 
9m), the gap function can be expanded as 



A(q) 



A + WM ( q'y - —q'x 



(a) 




qy 0.0 



FIG. 5: (a) Contours of the gap function Ak = [C+(k) — 
5-(k)]/2 for Pa = 5.4 kbar in CO (II) phase where < ky < 
2tt and < fc^; < 2n. The two Dirac points emerge from 
the M point kj\/ = (7r,±7r). Points O, X and Y denote (0,0), 
(tTjO) and (0,7r) , respectively, (b) The same contour Aq in the 
vicinity of the M point, which is taken as the origin (q = k — 
kA/). The red curve represents the energy contour Aq = A. 
The two red lines denote the angles 9min and Omax defined in 
the text l|4.18|l . Note that on Fig. (b), the scales are different 
along the two axes. 



A 



dct Sm ,2 
Qx 

WM 



WAiq'y 



(4.24) 



This explicitly shows that dctSM is the single parameter 
for the merging condition. 

The Pa-dependence of A, Vj, w''^ and det5'M for Va = 
0.18 eV and V^ = 0.05 eV are shown in Fig. H The 
merging pressm'e, calculated from detS'A/ = 0, is found 
to be P^ ~ 5.07 kbar, which almost coincides with that 
obtained from the direct calculation of energy bands. 
The merging axis g^-axis is very close to the qy axis, 
Om = 1-73 as shown in Fig. [5lja). 



D. Relation between effective Hamiltonians for 
a-(BEDT-TTF)2l3 and two-component systems 

In this section, wc compare the physics of the emer- 
gence of Dirac points in a-(BEDT-TTF)2l3, with a pos- 
sible similar transition in two-component systems such as 
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FIG. 6: Pressure dependence of (a) A, (b) Vj, (c) lOg', (d) Wg', 
(e) det^M and (f) dma,^, Ormr, for Va = 0.18 eV and 14 = 0.05 
eV. The pressure Pa is in kbar. 



graphene or boron nitride (BN) under stress. Near the 
M-point, the electronic spectrum of these two systems 
can be described by the following Hamiltonian^ii^iSi: 



^BN 



(q) 
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Afcra 



M 
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^*+ 9^ ^ *"^9^ 



-M 



2m 



vqx(j2 



(4.25) 



In graphene Af = 0, and in boron nitride A/ 7^ is 
proportional to the energy difference between sites B and 
N. The parameter A* drives the transition from a phase 
with two Dirac points (A* < 0) separated from 5q = 
2^— 2mA* to an insulating phase with a gap A* > 0. 
This transition has been studied in refs.— 1^— . 

Here we show that this transition has the same struc- 
ture as the transition described above in a-(BEDT- 
TTF)2l3. To do this, we have to compare the Luttingcr- 



Kohn Hamiltonian (j4.15p and the " boron-nitride" Hamil- 
tonian (|4.25p . This is possible after a rotation of angle 
—ip where 

tan^ = ^, (4.26) 

from which we obtain a new Hamiltonian: 

/iij^.(q)=exp'*"^/iBjv(q)exp"'*'^^ . (4.27) 

This Hamiltonian has the form 



'iBwlq) 



(A + w^ql)(7^ + Wiqldi + vqx(T2 



with the parameters uii, w^ and A given by 
1 1 M 

Wi = - — COSip 



2m 



2m ^A.p + A2 



1 . 1 A, 

W3 = - — sm Lp = 



2m 



(4.28) 



(4.29) 



(4.30) 



2m ^AP + A2 
A = v/a!Tm2. (4.31) 

This Hamiltonian therefore appears as a peculiar ease 
of the general Hamiltonian (|4.15p , with the following pa- 
ws and all other 
9, wf = wz sin^ 0, 



rameters A, 



,,vv 



vv 

Wl, Wg 



parameters being zero (w^ = wi sm" 

ve = VxCOsO). In this case, the determinant of the sta- 

2 

bility matrix is simply dot Sm = 2A^3' T^^^ merging 
transition has the same structure as the one described 
above. Here we have simply 
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In this simple case, the direction of the merging stays 
constant and perpendicular to the velocity (see Fig. [T]). 
The BN Hamiltonian is interesting since it explicitly re- 
veals the role of the different parameters. Whereas w^ 
drives the topological transition, the parameter wi con- 
trols the opening of a gap in both phases. The situation 
in Q;-(BEDT-TTF)2l3 is more involved since there are 9 
parameters but the main features of the topological tran- 
sition can be understood within a simplified Hamiltonian 
with only 4 parameters. 



qy 0.0 




FIG. 7: Contours of A(q) for the BN model. Here we have 
taken the same parameters as for the CO(I)-CO(II) transition 



in a(BEDT)-TTF)2l3, A, wi 



w-j = wf^, Vx, other w^ 



■ujg , Vy being zero. In this simple case, since the parameter 
uJg^ = 0, the Dirac points stay along the y-axis. The red 
curve corresponds to A(q) = A. The two lines indicate the 
directions O^ax and Omin- 



V. BERRY CURVATURE ASSOCIATED WITH 
THE TOPOLOGICAL TRANSITION 



In this section we show that the emergence of the Dirac 
pair is well characterized by the appearance of a non triv- 
ial Berry curvature i3i-|-(k) for the conduction band. This 
Berry curvature, which plays the role of a k dependent 
effective magnetic field in the Brillouin zone, appears to 
be sharply peaked around the two massive Dirac points 
(see Fig. [5]Jb)). We then show how to associate a Berry 
phase to each Dirac point. 



A. Full Brillouin zone computation of the Berry 
curvature within the four-band model 



For a multiband system, the Berry curvature i?„ (k) of 
the v}^ band is given hy^: 



B„(k) =rotkA„(k) 



(5.1) 



where A„ (k) denotes the so-called Berry connection and 
is written as 



A„(k) = -i < n(k)|ak|n(k) >, 



(5.2) 



and |n(k) > is an eigenvector of Eq. (|2.8p . 

For our model Hamiltonian (j2.5p which is diagonal in 
spin index and with non degenerate bands S,na{}^) that 
never cross for any k, the Berry curvature B„o-(k) = 
-B„o-(k)Uz can be computed from 

u n\ \^ <ma(kX„a(k)-fc.c. 



tn^n 



(Cn<.(k)-U.(k))2 ' 



where 



= Ea/3 dnaaOi)*dmPa(k)dk,_^eaflik) 



(5.4) 



are the nm intcrband matrix elements of the velocity 
operator. n,m = 1,2,3,4 are Blocli bands indices and 
a = A, A', B, C are sites in a-(BEDT-TTF)2l3. Note 



that since 



x,y 



r(k) 



\mna \ 



^(k)) , we immediately obtain 
that -B„<^(k) = -I]m5^„-Sm<T(k). Therefore computing 
_Bicr(k) fully characterizes our 3/4 filled system. 

Fig. |5] shows the typical result for i5i-|-(k) in the insu- 
lating CO(I) {Pa = 4.4 kbar) and CO(II) {Pa = 5.4 kbar) 
phases. In the CO (II) phase, there is a pair of massive 
Dirac particles and consequently the Berry curvature is 
sharply peaked around the position of the Dirac points 
k±. The peaks at kj- are strongly anisotropic and have 
the same magnitude but opposite sign, reflecting the so- 
called opposite chirality associated to the two Dirac par- 
ticles. Conversely, in the CO (I) phase, the intensity of 
the Berry curvature becomes small owing to cancelation 
of the positive and negative peaks. 



B. Berry curvature of the effective tw^o bands 
Luttinger-Kohn model 

In order to understand in a more quantitative way 
these numerical results, we show here how the Berry cur- 
vature Bi^(k) can be computed directly from the low 
energy effective 2x2 Luttinger-Kohn reduced Hamilto- 
nian iJi^^(kM + q) = /o(q)fo + f (q)'f. f is the vector of 
components (/i, /2, /a) (Eq. 14. 4p . For this 2x2 Hamilto- 
nian, the Berry curvature i?if (q = k — kM) (in cartesian 
and polar coordinates) reads 



^it(q) 
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(5.5) 



2QA3(q) ■ 
From Eq. (j4.15p we obtain: 

A(2wf -We' - v^wf')q + ve{wlwf' - wfw^')q^ 



Bn{q,e) 



2[{A + wlq^)^ + v^q^ + {wf)^q' 



i3/2 



where we have defined the derivatives 



(5.6) 



wf = {w\^ 
wf = {wf 
vg' = —Vx sin + Vy cos ( 



wf'') sm2e + 2w^y cos29, 

wf^ ) sin 26* -h 2^3*' cos 26*, (5.7) 




The expression Eq. (|5.6p shows that the wf F term is es- 
sential to have a non- vanishing Berry curvature. We note 
that this expression is quite general. In particular, for the 
boron nitride Hamiltonian described in previous section, 
we obtain 



S""(q) = 



-Awiv^qy 



((A 



wsq'y)^ 



vlql 



n^ni 



it)' 



(5.8) 



The q^ term in the numerator (that determines the large 
q tail of the Berry curvature) vanishes here because wf 
and lUg have identical 6 dependences. 



C. Berry phase associated to Dirac points 

The integral of the Berry curvature i?i-t-(k) over the 
full Brillouin zone (BZ) is a topological quantity called a 
Chern number— 



Ch 
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2^ 



BZ 



d5Si^(k). 
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FIG. 8: (a) Berry curvature -Bit(k) in CO(I) with Pa = 4.4 
kbar, (b) Berry curvature -Bi-f(k) in CO(II) with Pa = 5.4 
kbar, where Va = 0.18 eV and 14 = 0.05 eV. The center 
denotes the M point. 



In our system, Chi^ = since the curvatures associated 
to the two Dirac points exactly compensate. In order 
to characterize the contribution of each Dirac point, we 
define a Berry phase which is the integral over an appro- 
priate region in k space around a given Dirac point as 
explained in^l: 
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(5.10) 
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FIG. 9: Pa-dependence of the Berry phase ir(k±)[ defined 
in Eq. (|5.10[) (soUd fine) compared with the expression r± 
given in Eq. (IS.llfl (dashed fine). The inset denotes <?+„, the 
y-component of one Dirac point k+, where q = k — k*/. 



The area 5'(k-|-) denotes the region around a given Dirac 
point k-t where A(k) contours are closed (A(k) < A). 
This is the region inside the red curve in Figures [5] and 



m 



Here we compare r(k-i-) with the following expression 



r± 



^l/,_A(k^ 
2 V A(kM) 



(5.11) 



This expression was obtained in Ref. |2J for a general 
class of 2 X 2 Hamiltonians that include the BN model 
for which we find 



r.^4 1 



Wi 



^■^1+' 



4ii 



M 



2 \ ^AP + A2 



(5.12) 
Therefore in the vicinity of the transition (ws — > 0_), V± 
varies as 

|r,|4^. (5.13) 

At the same time, the separation 8q between Dirac points 
varies as 



5g: 



(-AW3) 



1/2 



(5.14) 



The Pa-dependence of r(k±) is compared with that of 
r± in Fig. m Although the behavior of r(k±) is qualita- 
tively similar to F-t, there is a slight difference between 
these two quantities, which may come from the effect of 
the other bands in Eq. (|5.3p . This multi-band effect is 
discussed in a separate paper—. Both quantities F(k±) 
and F-t vanish below a critical pressure which coincides 



with the merging pressure P^ found in Fig [6l within 
numerical accuracy. This is consistent with the effec- 
tive Hamiltonian which leads to |F(k±)| oc {Pa - P<f )^ 
and q^y oc [Pa — P^)^/'^. This is also seen quite simply 
in the BN model for which we obtain, from Eqs. (|5.13p 

and (|5.14p . |r±| ex w§ ex (det^A/)^ and g+y 
(det^M)^/'. 



,1/2 



VI. SUMMARY AND DISCUSSION 

In this paper, we re-examined the electronic proper- 
tics in the stripe charge ordered states in the organic 
conductor a-(BEDT-TTF)2l3. With increasing pressure, 
a transition occurs from an insulating charge ordered 
(CO) state to a metallic charge ordered (COM) state. 
In this work, we have found a new topological tran- 
sition which further splits each of the CO and COM 
phases into two phases, respectively CO(I)-CO(II) and 
COM(I)-COM(II). We specifically consider the CO(I)- 
CO(II) transition which separates a phase (CO (I)) with 
a usual gap at the M-point of the Brillouin zone from a 
phase (CO (II)) where the gap exhibits a double minimum 
with a local maximum at the M-point. This transition 
corresponds to the emergence of a pair of Dirac points. 

This modification of the band structure is described 
by an effective 2x2 low energy Luttinger-Kohn Hamil- 
tonian with 9 parameters that can be extracted from a 
numerical Hartree calculation. From a detailed study 
of this Hamiltonian and the corresponding energy spec- 
trum, we show that the transition is driven by a single 
quantity detS'M (Eq. (|4.12p ) which is an appropriate com- 
bination of these parameters and whose sign changes at 
the transition. A similar scenario occurs for the COM(I)- 
COM(II) transition. The Hartree contribution, which 
induces on-site potentials, is crucial to break inversion 
and time-reversal symmetries, leading to the CO state. 
The exchange (Fock) term modifies the hopping ener- 
gies but does not break the symmetry, so that it is not 
responsible for the CO state. It may only modify the rel- 
ative stability of the different phases. Moreover, in the 
ZGS phase, it does not disturb the Dirac spectrum as 
shown by first principle calculations^'^' ^'*. Similarly we 
believe that going beyond the mean field approximation 
does not change qualitatively the topological transition 
around the M point but may only modify the parameters 
of the Luttinger-Kohn Hamiltonian. 

We compare the structure of this transition with a sim- 
ilar situation which occurs in a simple model for boron- 
nitride^— i^, where the emergence of a pair of Dirac 
points is quite comparable, although the driving forces 
are different. 

The existence of a pair of Dirac points is characterized 
by a special structure of the the Berry curvature inside 
the Brillouin zone. In the CO (II) and COM(II) states, 
the latter shows two sharp peaks with opposite signs. 
On the other hand, in CO(I) and COM(I) phases, the 
Berry curvature becomes very small owing to cancelation 



11 



of the positive and negative contributions. The existence 
of the Dirac point is also verified by integrating the Berry 
curvature over a region hmited by a closed energy contour 
around a single point. 

It would be interesting to observe directly the modifica- 
tion of the energy spectrum associated to the topological 
transition by angle-resolved photoclcctron spectroscopy 
experiments. Moreover this topological transition could 
be probed in a magnetic field by the modification of the 
Landau level structure^, therefore by e.g. magnetoresis- 
tance and Hall experiments, in a range of pressure which 
is accessible experimentally. Finally we mention the ex- 
istence of another conductor, a-(BETS)2l3 which is also 
a good candidate for a Dirac-like electron spectruni^. 
Moreover, since the BETS molecule contains relatively 
heavy atoms, the spin-orbit interaction may be effective 
leading to non-zero Chern number as shown in the Weyl 
Hamiltonian with the spin-orbit interactioi*^. 
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pair of minima of the gap in the f spin band as shown in 
Fig. [TU] The gap-minima of the f spin band are located 
at kO| ^ ±(0.957r,-0.7l7r) for Pa = 5.4kbar, K = 0.18 
eV and H = 0.05 eV. 

(3) ky/7i 

0.0 



kx/ii 




kx/n 



Appendix A: Band Structure in the stripe Charge 
Ordered state 

In the stripe CO state, the energy band with f spin is 
different from that with ^ spin owing to inversion symme- 
try breaking and time reversal symmetry breaking. For 
the CO(I) phase ( Va = 0.17 eV, and Pa = 5 kbar in 
Fig. [2pS, the minimum point of the energy gap is lo- 
cated at kO-|- = kM = (tt, ±7r) (the M-point ) for ^it(k), 
and at kOj, = ±(0.367r,-0.627r) for Cu(k). The t spin 
band has a small gap A-|- = 0.0033 eV at the M-point 
and is very anisotropic in its vicinity. The J, spin band 
has a large gap Aj^ = 0.0171 eV and is more isotropic. In 
the CO (II) phase (Fig. [2]), on the other hand, there is a 



FIG. 10: The energy bands of ^lo- and ^2ct for '|" spin band 
(a) and the 4, spin band (b) in the CO(II) for Pa = 5.4kbar, 
Va = 0.18 eV and Vb = 0.05 eV. The center denotes tlie M 
point. 



The energy dispersion in COM (II) is similar to that of 
CO(II). However, for COM(II), the Fermi level is located 
in the conduction band for the f spin band, and in the 
valence band for the | spin band (Fig [3] d) . When the 
pressure is decreased, those two minimum points of the 
t spin bands merge at the M-point while the quantitiy 
Cit ~ ^2t ( ?i^ ~ ^2^ ) remains finite on the boundary line 
between COM(I) and COM(ll). 



N. Tajima and K. Kajita K, Sci. Tech. Adv. Mater. 
024308 (2009). 



10, 



H. Seo, C. Hotta, and H. Fukuyama H, Chem. Rev. 104, 



12 



5005 (2004). 

A. Kobayashi S. Katayama, and Y. Suzumura, Sci. Tech. 

Adv. Mater. 10, 024309 (2009). 

* H. Kino and H. Fukuyama, J. Phys. Soc. Jpn. 64, 4523 
(1995). 

^ H. Seo, J. Phys. Soc. Jpn. 69, 805 (2000). 
^ C. Hotta, J. Phys. Soc. Jpn. 72, 840 (2003). 
■^ T. Takahashi, Synth. Met. 133-134, 26 (2003). 

* N. Tajima, A. Ebina-Tajima, M. Tamura, Y. Nishio, K. 
Kajita, J. Phys. Soc. Jpn. 71, 1832 (2002). 

® A. Kobayashi, S. Katayama, and Y. Suzumura, J. Phys. 

Soc. Jpn. 74, 2897 (2005). 
^° K. Kajita, T. Ojiro, H. Fujii, N. Nishio, H. Kobayashi, A. 

Kobayashi, R. Kato, J. Phys. Soc. Jpn. 61, 23 (1992). 
^^ S. Katayama, A. Kobayashi, and Y. Suzumura, J. Phys. 

Soc. Jpn. 75, 054705 (2006). 
^^ R. Kondo S. Kagoshima, and J. Harada, Rev. Sci. Instrum. 

76, 093902 (2005). 
^^ H. Kino and T. Miyazaki, J. Phys. Soc. Jpn. 75, 034704 

(2006). 
^"^ S. Ishibashi, T. Tamura, M. Kohyama, and K. Terakura, 

J. Phys. Soc. Jpn. 75, 015005 (2006). 
^^ A. Kobayashi, S. Katayama, Y. Suzumura, and H. 

Fukuyama, J. Phys. Soc. Jpn. 76, 034711 (2007). 
^^ M.O. Goerbig, J.-N. Fuchs, G. Montambaux, and F. 

Piechon, Phys. Rev. B 78 , 045415 (2008). 
^^ A. Kobayashi, Y. Suzumura, and H. Fukuyama, J. Phys. 

Soc. Jpn. 77, 064718 (2008). 
^^ N. Tajima, S. Sugawara, R. Kato, Y. Nishio, and K. Kajita, 

Phys.' Rev. Lett. 102, 176403 (2009). 
^^ K. Morinari, T. Himura, and T. Tohyama, J. Phys. Soc. 

Jpn. 78, 023704 (2009). 
^° P. Dietl, F. Piechon, and Montambaux G Phys. Rev. Lett. 

100, 236405 (2008). 



^^ G. Montambaux, F. Piechon, J.-N. Fuchs, and M.O. Go- 
erbig, Eur. Phys. J. B 72, 509 (2009). 

^^ G. Montambaux, F. Piechon, J.-N. Fuchs, and M.O. Go- 
erbig, Phys. Rev. B 80, 153412 (2009). 

^^ M. V. Berry, Proc. Roy. Soc. London A 392, 45 (1984). 

^^ J.-N. Fuchs, F. Piechon, M.O. Goerbig, and G. Montam- 
baux, Eur. Phys. J. B 77, 351 (2010). 

^^ T. Mori, A. Kobayashi, Y. Sasaki, H. Kobayashi, G. Saito, 
and H.Inokuchi, Chem. Lett. 957 (1984). 

2^ T. Mori, H. Mori, and S. Tanaka, BuU. Chem. Soc. Jpn. 
72, 179 (1999). 

^'^ T. Kakiuchi, Y. Wakabayashi, H. Sawa, T. Takahashi, and 
T. Nakamura, J. Phys. Soc. Jpn. 76, 113702 (2007). 

^^ Y. Takano, K. Hiraki, T. Takahashi, and H.M. Yamamoto, 
J. Phys. Soc. Jpn. 79, 104704 (2010). 

■^^ A. Kobayashi, S. Katayama, K. Noguchi, and Y. Suzu- 
mura, J. Phys. Soc. Jpn. 73, 3135 (2004). 

^° A. Kobayashi, S. Komaba, S. Katayama, and Y. Suzu- 
mura, J. Phys., Conf. Ser. 132, 012002 (2008). 

^^ J.M. Luttinger and W. Kohn W Phys. Rev. 97, 869 (1955). 

^^ S. Katayama, A. Kobayashi, and Y. Suzumura, Eur. Phys. 
J. B 67, 139 (2009). 

^^ The possibility of a transition towards a local maximum 
is also possible, but does not correspond to the physical 
situation discussed in this paper. 

'^'* D.J. Thouless, Topological Quantum Numbers m Nonrela- 
tivistic Physics (World Scientific, Singapore, 1998) (1998). 

^^ Y. Suzumura and A. Kobayashi, private communication. 

^•^ K. Hiraki, S. Harada, K. Arai, Y. Takano, T. Takahashi: 
J. Phys. Soc. Jpn. 80, 014715 (2011). 

^'^ C.L. Kane and E.J. Mclc, Phys. Rev. Lett. 95, 226801 
(2005). 



